LongCART (version 3.2)

StabCont : parameter stability test for continuous partitioning variable

Description

Performs parameter stability test (Kundu and Harezlak, 2019) with continuous partitioning variable to determine whether the parameters of linear mixed effects model remains same across all distinct values of given continuous partitioning variable.

Usage

StabCont(data, patid, fixed, splitvar)

Arguments

data

name of the dataset. It must contain variable specified for patid (indicating subject id) and all the variables specified in the formula and the StabCont(data, fixed, splitvar)partitioning variable of interest specified in splitvar.

patid

name of the subject id variable.

fixed

a two-sided linear formula object describing the fixed-effects part of the model, with the response on the left of a ~ operator and the terms, separated by + operators, on the right. Model with -1 to the end of right side indicates no intercept. For model with no fixed effect beyond intercept, please specify only 1 right to the ~ operator.

splitvar

the continuous partitioning variable of interest. It's value should not change over time.

Value

p

It returns the p-value for parameter instability test

Details

The continuous partitioning variable of interest. It's value should not change over time.

$$Y_i(t)=W_i(t)theta+b_i+epsilo{n}_{it}$$

where $W_i(t)$ is the design matrix, theta is the parameter associated with $W_i(t)$ and b_i is the random intercept. Also, $epsilo{n}_{it}N(0,sigm{a}^2)$ and $b_{i}N(0,sigm{a}_u^2)$. Let $X$ be the baseline continuous partitioning variable of interest. StabCont() performs the following omnibus test

$H_0:thet{a}_{(g)}=thet{a}_0$ vs. $H_1:thet{a}_{(g)}^{=}thet{a}_0$, for all g

where, $thet{a}_{(g)}$ is the true value of $theta$ for subjects with $X=C_g$ where $C_g$ is the any value realized by $X$.

References

Kundu, M. G., and Harezlak, J. (2019). Regression trees for longitudinal data with baseline covariates. Biostatistics & Epidemiology, 3(1):1-22.

See Also

StabCont, LongCART, plot, text

Examples

# NOT RUN {
#--- Get the data
data(ACTG175)
               
#--- Run StabCont()                
out<- StabCont(data=ACTG175, patid="pidnum", fixed=cd4~time, splitvar="age")
out$pval
# }